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Abstract 

We have constructed the surface phase diagram of the BaTi 03 ( 001 ) surface by employing an evo¬ 
lutionary algorithm for surface structure prediction, where the ferroelectric polarization is included 
as a degree of freedom. Among over 1000 candidate structures explored, a surface reconstruction 
of (2xl)-TiO is discovered to be thermodynamically stable and have the p2mm plane group sym¬ 
metry as observed experimentally. We find that the influence of ferroelectric polarization on the 
surface free energy can be either negligibly small or sizably large (over 1 eV per (2 x 1) supercell), 
depending strongly on the surface structure and resulting in a significant distinction of surface 
phase diagram with varying ferroelectric polarization. It is therefore feasible to control the surface 
stability by applying an external electric field. Our results may have important implications in 
understanding the surface reconstruction of ferroelectric materials and tuning surface properties. 

PACS numbers: 68.35.B-, 77.80.-e, 05.70.-a, 71.15.Mb 
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The search for stable surface structures is a key subject of surface science and of great 
importance to fundamental research as well as practical applications, like photovoltaics, 
catalysis, and sensors jl-d]. Density functional theory (DFT) in combination with ab initio 


Kl. 


In 


thermodynamics is an indispensable tool because of the atomic insight it provides 
this approach, the surface free energy is expressed as a function of stoichiometry and atomic 
chemical potentials so as to consider the varying growth conditions, and the minimization 
of the surface free energy predicts stable surface phases. Very recently, the approach has 
been developed for semiconductors to include the electron chemical potential as a new pa¬ 
rameter, which can be generally applied to study the effects of bulk dopants on properties 
of semiconductor surfaces and interfaces j^. Further generalization of the approach to other 
systems would be interesting. 

In ferroelectric materials, the ferroelectric polarization couples strongly with the crystal 
structure, and consequently any change of the ferroelectric polarization will in turn affect 
the structural stability. Thus, the ferroelectric polarization that can be easily controlled by 
external electric field is an important degree of freedom for ferroelectric surfaces. However, 
ferroelectric polarization has hardly been taken into account in the previous calculations 
of surface phase diagrams of ferroelectric materials. Recent first-principles study indeed 
showed that surface stability of ferroelectric lithium niobate is different for the positively 
and negatively polarized surfaces, which is actually driven more by the different surface 
termination than intrinsic ferroelectric polarization [8]. 

As a prototypical ferroelectric material, barium titanate (BaTiOs, BTO) plays a vital 
role in numerous applications and has been intensively studied theoretically { q - 16] As shown 
in Fig. 1, BTO has a perovskite structure and undergoes a structural transition from high- 
symmetry cubic phase to low-symmetry tetragonal phase when lowering temperature across 


-400 K 


0 . 


The (001) face is a stable cleavage plane and has rich surface reconstructions, 
including (1x1), (2x1), c(2x2),(2x2),(v^x \/5), (3x1), (3x2), and (6x1) periodicities lisl- 
23]. Among them, the (2 x 1) reconstruction has recently attracted much attention jl^. bJ. 
251 ]. Two different surface structure models have been proposed. However, one structure 
model [^, 24] displays the pm plane group symmetry, in contradiction with the w2mm 
symmetry identified by recent x-ray diffraction experiments [2^; the other one [2^ has 
the correct symmetry but is energetically less stable than the former. The actual atomic 
structure of the BTO (001)-(2 x 1) surface remains elusive. 
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FIG. 1. (Color online) The atomic strnctnre of bnlk BaTiOs in paraelectric and ferroelectric phase. 
The polarizations and are determined by the displacement between O and Ba/Ti along [001]. 
The green, blue and red balls represent Ba, Ti and O atoms, respectively. 


In this work, we will consider ferroelectric polarization as an extra degree of freedom 
to calculate surface phase diagram of ferroelectric materials. Specifically, we perform first- 
principles calculations for the BTO(OOl) surface, focusing on (2 x 1) as well as (1 x 1) 
reconstructions to find stable surface configurations and to reveal the effects of ferroelec¬ 


tric polarization. By employing a newly designed evolutionary algorithm [2fi| for efficiently 
searching (meta)stable configurations and calculating over 1000 possible structure models, 
we predict a surface phase diagram containing many new surface structures, including a 
thermodynamically stable (2 x l)-TiO phase that has the p2mm plane group symmetry 
observed experimentally (251]. More importantly, we find that the influence of ferroelectric 
polarization on the surface free energy can be either negligibly small or sizably large (over 
1 eV per (2 x 1) supercell for BTO(OOl)), depending strongly on the surface structure. As 
a result, the surface phase diagram changes significantly with varying ferroelectric polariza¬ 
tion. These findings suggest a unique way to control surface structures and properties of 
ferroelectrics. 

First-principles calculations were performed with DFT as implemented in the Vienna 


ab initio simulation package 


27l |. using the projector augmented wave methodj21 


and the 


Perdew-Burke-Ernzerhof (PBE) exchange correlation functional |22l, l2Q|- The method for 
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predicting surface reconstructions was based on an evolutionary algorithm as implemented 
in the USPEX package j^, allowing variable-composition structure searches, where the num¬ 
bers of atoms in the surface region are varied to yield the global minimum of the surface free 


energy (see Refs. 


26, 


32, 


for more details). To demonstrate the effect of ferroelectric po¬ 


larization on surface reconstructions, we focused on discussing the Ti02-terminated surfaces, 


which have been extensively observed experimentally and studied theoretically 






and considered only an ideal bulk-terminated phase for the BaO-terminated surface. The 
BTO(OOl) surfaces were modeled by periodic slabs composed of four -Ti02-Ba0- bilayers 
plus a Ti02 termination together with a 15 A thick vacuum layer [3^. Different surface 
stoichiometries were considered by adding a layer of Ti^Oy {x=0, 1, 2, y=0, 1, 2, 3, 4) in 
a (2 X 1) surface supercell on the otherwise ideal Ti02-terminated surface. The bottom 
three bilayers were fixed at their bulk configuration, and the other layers were relaxed us¬ 
ing the conjugate gradient algorithm till residual forces were smaller than 0.01 eV/A^. A 
Monkhorst-Pack A:-point mesh with reciprocal-space resolution of 2 tt x 0.03 A“^ and a 400 
eV planewave cutoff energy were used. Dipole correction was employed in slab calculations 
for removing artificial interactions between the slab and its periodic images. 

Thermodynamical stability of a surface structure is determined by the surface free energy, 
7 = G'siab - Gref “ AnBa/^Ba “ AuTi/iTi “ AnoyUo, where Gsiab and Gref are the Gibbs free 
energies of the slab and the reference system that was selected as the ideal Ti02-terminated 
surface. Anea, Auxi and Auq denote the changes in the number of atoms with respect to the 
reference system. All these quantities correspond to a (2 x 1) surface supercell if not specified 
otherwise. pBa, hTi and po are the atomic chemical potentials. The accessible boundary of 
chemical potentials is defined by thermal equilibria between bulk BaTiOs and other phases, 
including bulk Ba, bulk Ti, bulk BaO and bulk Ti02. Herein we approximated the Gibbs free 
energy by the DFT total energy, excluding the vibrational contribution. The approximation 
has been found to be satisfactory for our study: the phase diagram qualitatively remains 
unchanged when the temperature effect was considered (see the Supplemental Material |35|). 

At first we exclude the contribution of ferroelectric polarization by fixing the lower three 
-Ti02-BaO- bilayers at the cubic bulk structure (i.e., paraelectric phase), as typically done 


in previous studies 


Q,y, 


25|. In contrast to previous studies, we coinputed more sur¬ 


face configurations (over 1000) using an advanced evolutionary algorithm |2fil|. and obtained 
many new stable surface structures. Figure 2(a) shows the computed surface phase dia- 
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FIG. 2. (Color online) (a) Surface phase diagram of paraelectric BTO(OOl) for (1 x 1) and (2 x 1) 
reconstructions without the effect of ferroelectric polarization. The red solid lines bound the 
chemical stability ranges of BTO. The precipitation lines of bulk BaO, Ti and Ti02 are labeled as 
1, 2 and 3, respectively, which bound the accessible chemical potential range defined by thermal 
equilibria, (b) Schematic atomic structures of stable surface phases. The blue and red balls 
represent Ti and O atoms, respectively. 


gram of BTO(OOl) for (1 x 1) and (2 x 1) reconstructions, and Fig. 2(b) shows the atomic 
configurations of the stable phases. It can be seen from Fig. 2(a) that the (1 x 1) ideal BaO- 
terminated surface is stable at 0-rich and Ba-rich conditions. With decreasing /tq and //Ba, 
other phases become increasingly more stable. These stable phases include a double-layer 
Ti02-termination model (2 x 1 )-Ti 204 [3, 24], a TiO adunit model (2 x 1)-Ti0 formed by 
adding a TiO unit vertically at the hollow site, and two Ti adatom models, (2 x 1)-Ti and 
(1 x 1)-Ti, formed by adding a Ti atom at the hollow site in the surface supercell. Note 
that another double-layer model (1 x 1)-Ti02 in Fig. 2(b) (not shown in the phase diagram) 
has surface free energy very close (within ~ 2 meV/A^) to that of the (2 x 1 )-Ti 204 phase. 
Similar geometrical features of Ti=0 titanyl are found in these two thermodynamically 
degenerate phases. 

In experiment, the BTO(OOl) (2 x 1) surface reconstruction is obtained by Ar"*" ion 
sputtering and subsequently annealing native (1x1) surface |25|. Our surface phase diagram 
suggests three (2 x 1) phases, including the known double-layer model (2 x 1 )-Ti 204 |18l.l24|. 


Another double-layer (2x1) model proposed by previous work [2^ does not appear in the 
phase diagram for the reason that it is thermodynamically less stable and its surface free 
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energy is far higher (~ 1.4 eV) than that of the donble-layer model (2 x 1 )-Ti 204 shown 
in the phase diagram. Carefnl analysis of X-ray diffraction data indicates that the (2 x 1) 
surface has the p2mm plane group symmetry |25|. Based on this information, the (2 x 1)- 
Ti 204 and (2 X 1)-Ti models that have the pm symmetry are excluded. The only remaining 
(2 X 1) model with the p2mm symmetry is (2 x 1)-Ti0, which is a thermodynamically stable 
phase located within the chemical stability ranges of BTO [see Fig. 2(a)]. Interestingly, 
among over 1000 structures suggested by the evolutionary algorithm, the (2 x 1)-Ti0 model 
is the only one that satisfies both conditions of energy and symmetry. We thus attribute the 
formation of (2 x 1) reconstruction to the addition of TiO units, which is consistent with 
recent experimental observation of TiO adunits on the c(2x2) BTO(OOl) surface (36|. 


We then consider different ferroelectric polarizations {Pi and Pi) by fixing the lower 
three -Ti02-Ba0- bilayers at the tetragonal bulk structures (ferroelectric phase) as shown 
in the right panel of Fig. 1. Figure 3 shows the calculated surface phase diagram for two 
opposite ferroelectric polarization orientations. The interesting (2xl)-TiO structure remains 
as a stable phase in the surface phase diagram for both types of polarizations. However, 
distinct variations of phase diagrams in different polarization conditions can be found, e.g., 
(2xl)-TiO is stable in the Pi condition over a remarkably extended chemical potential 
range, while under the condition (2xl)-TiO becomes unstable unless the po is fairly 
low. The distinction between the phase diagrams under different polarizations indicates 
that at a certain experimental circumstance (e.g., the oxygen chemical potential is within 
—3 eV ~ —5 eV), the thermodynamically stable phases are different (i.e., (2xl)-TiO and 
(1x1) ideal BaO-terminated surfaces). Thus the external electric field induced ferroelectric 
switching of the substrate might lead to appearance of different surface structures if kinetic 
factors were not considered. 


To understand this phenomenon, we quantitatively analyze the influence of ferroelectric 
polarization on the phase diagram by calculating the relative surface Gibbs free energies, 
A 7 (P 4 ^/-|-)= 7 (P 4 ^/-|-)—7(P=o), for various surface structures under two opposite polarization 
conditions (the case without ferroelectric polarization is taken as the reference). As shown 
in Table 1, the obtained nonstoichiometric reconstruction phases can be divided into two 
types: the ones with a Ti 02 overlayer and the other ones with a Ti^;©^; {y<2x) adunit on the 
primary Ti02 termination. The surfaces of the “adunit”-type show a considerable energy 
difference of /2^^{Pi) — A 7 (P|) ( ~ 1.0 eV/(2xl)cell), indicating significant influence of the 
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FIG. 3. (Color online) Same as Fig. 2(a) except that the effect of ferroelectric polarization is 
included. Different ferroelectric polarizations, (a) Pi and (b) as illustrated in Fig. 1(b), result 
in significantly different surface phase diagrams. 


ferroelectric polarization on the surface stability. In contrast, the corresponding influence is 
much smaller for the “overlayer”-type surfaces. Our results of structural relaxation explicitly 
show that the detailed surface atomic structures of stable phases hardly change with the 
effect of ferroelectric polarization. Thus, the remarkable difference between A 7 (Pj^) and 
A'y{Pi) does not come from structural relaxation but is mainly caused by electrostatic 
interactions. 


Generally, termination of the spontaneous polarization of ferroelectric materials always 
gives rise to discontinuity of polarization at surfaces, leading to surface polarization charges 
(and surface metallicity) whose signs depend on the direction of polarization (see more details 


in Supplementary Material 


35]). As shown in Figs. 4(a) and 4(b), these surface polariza¬ 


tion charges generate a depolarization held (or an internal electric held), whose direction 
is opposite to that of the polarization. With a constant non-zero depolarization held, the 
electrostatic energy of ferroelectric surfaces would diverge with increasing thickness. Such 
an electrostatic instability, however, can be eliminated by compensating the depolariza¬ 




tion held through various mechanisms, like introducing an external electric voltage [; 
surface adsorption or surface reconstruction. For the intrinsic mechanism of surface recon¬ 
struction, the compensation of the depolarization held depends significantly on the direction 
of polarization, leading to distinct surface free energies A'y{Pi) and A'y(Pi). 










TABLE 1. Relative surface free energies [A 7 ] and Bader charges for the snrface adnnits on the 
Ti 02 termination for the non-stoichiometric phases in different polarization conditions. The case 
without ferroelectric polarization is taken as the reference of A 7 . 


Phase 

A 7 (eV) 


Charge (e) 

Pi 

Pt 

Pi 

Pt 

( 2 xl)-Ti 204 

0.47 

0.02 

0.03 

0.02 

(Ixl)-Ti02 

0.38 

-0.09 

0.06 

- 0.02 

(2xl)-TiO 

-0.82 

0.32 

1.18 

1.15 

(2xl)-Ti 

-0.87 

0.37 

1.46 

1.45 

(Ixl)-Ti 

- 1.00 

0.48 

2.58 

2.59 

( 2 xl)-Ti 20 

-0.77 

0.17 

1.50 

1.53 


Different surface reconstructions result in different surface electrostatic potentials and 
thus correspond to varying surface dipoles. For a specified surface reconstruction, if the 
depolarization field is compensated by the surface dipole [see Figs. 4(c) and 4(d)], the 
electrostatic energy would get lowered by this compensation and the whole system thus gets 
stabilized. In contrast, reversing the depolarization field would yield higher electrostatic 
energy and larger surface free energy. This could well explain the polarization dependent 
behaviors, as demonstrated below. While the magnitude of surface dipole is not easy to 
quantify, a qualitative analysis is possible for the present system, considering that the normal 
oxidation states of Ti and O are +4 and — 2 , respectively. 

For the so-called “adunit”-type surfaces, the adunit that has Ti/0 ratio larger than 1/2 
is chemically unsaturated. When the adunit binds with the substrate, electron transfer from 
the adunit to the substrate occurs, resulting in a positively charged adunit as shown in Fig. 
4(c). For the Pi condition, the charge transfer decreases the surface polarization charge 
and the charge-transfer induced dipole compensates the depolarization field, resulting in 
a negative Aj{Pi). We denote such a charge-transfer induced compensation as the ionic 


surface compensation mechanism, as used in previous work 


39|. The effect gets inverted for 


the Pf condition, leading to a positive Aj{Pi). All these features are consistent with the 
calculation data (see Table 1). To analyze the results in more details, we present the calcu¬ 
lated Bader charges in Table 1. Due to the lower coordination number of surface atoms, the 
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FIG. 4. (Color online) (a) / (b) The discontinuity of ferroelectric polarization (P 4 ,/Pf) at the surface 
induces surface polarization charges, resulting in a depolarization field opposite to the direction of 
polarization. (c)/(d) Schematic of ionic surface compensation mechanism. Charge transfer between 
the adsorbate and substrate induces a surface dipole, which compensates the depolarization field 
and thus lowers the electrostatic energy of the system, leading to lower surface free energy. 


calculated charge of the adunits, e.g., TiO, which have formed Ti=0 double-bonded titanyl 
groups, is about -|-1.18e. A large surface free energy difference of “adnnit”-type phases [> 
1.0 eV/(2xl)cell] is induced through the ionic surface compensation mechanism described 
above. The detailed spatial distribntion of charge transfer is given in the Snpplementary 
Material 0 - 

For the so-called “overlayer”-type phases in which the overlayer itself is chemically satu¬ 
rated, there exists tiny charge transfer between the overlayer and the substrate (see Table 
1), suggesting that the effect of ionic surface charge compensation is negligible. The atomic 
rnmpling of Ti 02 overlayer, with the O atoms at the surface all above the Ti atoms, con¬ 
tributes a downward snrface dipole Ps{-1)- Snch a kind of structnral rnmpling has also been 
predicted for the bare snrfaces of perovskites which can lead to a relatively low catalytic 


activity of the surface0]. The anti-parallel/parallel configuration of the Pg to the P^ and 
Pi of the substrate leads to the stabilization/destabilization scenarios, respectively. Thus, 


10 










the polarization to the surface is suppressed for the condition but enhanced for condi¬ 
tion, resulting in the higher A 7 . The resulting changes of electrostatic potential alignment 
and the surface electronic structure by different types of reconstruction are given in the 
Supplementary Materialj35|. The obtained (2xl)-TiO surface shows n type metallicity in 
both polarization conditions. 

It should be noticed that the contributions of charge transfer and structural rumpling are 
strongly entangled and cannot be clearly distinguished by direct calculations, especially in 
the “adunit”-type phases like (2xl)-TiO. Nevertheless, compared to the relatively smaller 
[A 7 (P 4 ,) — A 7 (P|)] of (Ixl)-Ti 02 case 0.4 eV) which is basically contributed by the 
structural rumpling, the considerably larger [ A 7 (P 4 ,) — A 7 (P|)] of (2xl)-TiO case 
—1.1 eV) can be mainly attributed to charge transfer (i.e., the ionic surface compensation 
mechanism). 

The above results as well as related physical mechanism clearly indicate that the ferro¬ 
electric polarization plays a significant role in the surface stability of ferroelectric materials. 
This reveals a new degrees of freedom to affect the growth of the surface: in addition to 
tuning the growth condition, e.g., substrate temperature and partial pressure of source, it is 
convenient and feasible to control the surface stability by applying an external electric field. 

In summary, we have constructed the surface phase diagram of ( 2 x 1 ) and ( 1 x 1 ) 
BTO(OOl) reconstructions by employing a surface structure prediction method based on 
evolutionary algorithm and exploring over 1000 candidate structures. We predict a surface 
phase diagram containing many new surface structures, including a thermodynamically 
stable (2xl)-TiO phase that has the p2mm plane group symmetry as observed experimen¬ 
tally. Critically, the ferroelectric polarization has been included as a new parameter of 
surface structure prediction. We find the surface phase diagram changes significantly with 
varying ferroelectric polarization due to the ionic surface compensation mechanism. The 
distinguishing feature of ferroelectrics is the polarization switching upon applying external 
electric field, thus the control over surface stability is feasible by applying electric field. The 
underlying physical mechanism is expected to be quite general. Our results may help in 
tuning surface structures and properties of ferroelectric materials. 
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